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T2K has performed the first measurement of inclusive charged current interactions on carbon at 
neutrino energies of --^1 GeV where the measurement is reported as a fiux-averaged double differential 
cross section in muon momentum and angle. The flux is predicted by the beam Monte Carlo and 
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external data, including the results from the NA61/SHINE experiment. The data used for this 
measurement were taken in 2010 and 2011, with a total of 10.8 x 10^^ protons-on-target. The 
analysis is performed on 4485 inclusive charged current interaction candidates selected in the most 
upstream fine-grained scintillator detector of the near detector. The flux-averaged total cross section 
is (crcc),^. = {6M ±0.l3{stat) ±0.84:{syst)) x 10"^^ nucTeln ^ ^^^^ neutrino energy of 0.85 GeV. 



I. INTRODUCTION 

The T2K (Tokai-to-Kamioka) experiment is a long 
baseline neutrino oscillation experiment [l[ whose pri- 
mary goals are to make precise measurements of the ap- 
pearance of electron neutrinos and the disappearance of 
muon neutrinos at a distance where the oscillation is max- 
imal for the neutrino beam energy. This analysis is moti- 
vated, in part, by the fact that improved precision in the 
oscillation analyses requires better knowledge of neutrino 
interaction cross sections. 

Charged current (CC) neutrino-nucleon interactions at 
neutrino energies around 1 GeV have been studied in the 
past predominantly on deuterium targets 0, Many 
modern neutrino oscillation experiments like SciBooNE, 
MiniBooNE and T2K use heavier targets rich in carbon 
(and/or oxygen). Nuclear effects are important for those 
targets and, consequently, the neutrino interaction cross 
sections are not well known. It is therefore important to 
measure and understand these interactions to minimize 
systematic uncertainties for the oscillation measurement. 

For this purpose, we present a flux-averaged double 
differential inclusive muon neutrino CC cross section on 
carbon as a function of muon momentum and angle. An 
unfolding method that corrects for detector resolution, 
efficiencies and backgrounds is used to extract the flux- 
averaged total cross section and differential cross section 
in muon momentum and angle. The cross-section extrac- 
tion technique requires a detailed Monte Carlo (MC) sim- 
ulation including the neutrino flux prediction fScc. Ill B|) . 
a complete neutrino interaction model (Sec. IIIip . full de- 
scription of the near detector and accurate simulation of 
the readout electronics. The neutrino event generator 
NEUT [3| is used in the calculations of the selection effi- 
ciency for signal events and the rate of background events 
in the analysis. Additionally, a different event generator, 
GENIE , is used for cross-checks and for fake data stud- 
ies. Finally, the results of this analysis are compared with 
the predictions of the NEUT and GENIE generators. 

In this paper, we first summarize the T2K experiment 
in Sec. [ill where we present a description of the off-axis 
near detector (Sec. Ill C|) and, in particular, the tracker 
region where the neutrino interaction target is located. 
Sec. mil is devoted to a detailed description of how neu- 
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trino interactions are simulated in T2K. Sec. |TV] summa- 
rizes the reconstruction tools and the event selection in 
this analysis. The method used to extract the cross sec- 
tion is explained in Sec. |V] and the results are presented 
in Sec. IVII The conclusions are given in Sec. I VIII 

II. T2K EXPERIMENT 

T2K consists of an accelerator-generated neutrino 
beamline, a near detector complex 280 m downstream 
of the neutrino beam target and a far detector, Super- 
Kamiokande (SK), located 295 km away at an angle of 
2.5 degrees from the axis of the neutrino beam. Neutri- 
nos are generated from the 30 GeV proton beam of the 
Japan Proton Accelerator Research Complex (J-PARC) 
located in Tokai-mura on the east coast of Japan. The 
near detector complex is composed of a detector on the 
axis of the neutrino beam, called INGRID, and a detec- 
tor located 2.5 degrees off-axis, in line with SK, called 
ND280. INGRID is used primarily to measure the beam 
profile and stability, and the ND280 detector is used to 
measure neutrino fluxes and neutrino interaction cross- 
section properties, as the inclusive CC cross section pre- 
sented in this paper. 

In the following three subsections, the experiment is 
discussed in detail. We provide an overview of the neu- 
trino beamline in Sec. Ill Al Sec. IIIBI describes the T2K 
flux simulation included in the MCs used to extract the 
cross section, which represents a major systematic uncer- 
tainty. Sec. Ill CI describes the off-axis detector, empha- 
sizing the tracking detectors and target composition, as 
these are essential for the cross-section calculation. We 
refer to [l| for a complete description of the T2K experi- 
ment. 



A. Neutrino beam 

Fig. [T] depicts the neutrino beamline and the near de- 
tectors. The neutrino beam is produced by protons ac- 
celerated to 30 GeV kinetic energy. The proton beam 
has 8 bunches (6 before June 2010) with a 581 ns spac- 
ing. The protons in the spill are extracted and directed 
toward a 91.4 cm long graphite target aligned 2.5° off- 
axis angle from Kamioka. The target is installed in- 
side a magnetic horn that collects and focuses the pos- 
itively charged mesons (mainly pions and kaons) gener- 
ated by proton interactions in the target. Two additional 
magnetic horns are used to further focus the charged 



Target & horns 



uc^a^ cunnei 




beamdump ii^m^^.^^^ 




> Super-K 

ND280 



INGRID-280 m off-axis 
on-axis 

FIG. 1. A schematic view of the T2K neutrino beamhne and near detectors. 



mesons before they enter a 96 m long steel decay volume 
filled viritli helium. The mesons decay predominately into 
highly boosted muons and muon neutrinos, which prop- 
agate roughly in the direction of the decaying mesons. A 
beam dump stops most of the particles in the beam that 
are not neutrinos. Some high-energy muons pass through 
the beam dump and are observed by the muon monitor, 
providing information used to track the beam direction 
and stability. The analysis presented in this paper uses 
the data taken before March 2011, comprising a total of 
10.8 X IQis protons-on-targct (POT). 



B. Neutrino flux prediction 

A detailed description of the neutrino flux predic- 
tion can be found in [SL A FLUKA2008 0, [1 and 
GEANT3.21/GCALOR [i, ^ based simulation mod- 
els the physical processes involved in the neutrino pro- 
duction, from the interaction of primary beam protons 
in the T2K target, to the decay of hadrons and muons 
that produce neutrinos. The simulation uses T2K proton 
beam monitor measurements as inputs. The modeling 
of hadronic interactions is re-weighted using thin target 
hadron production data, including recent charged pion 
and kaon measurements from the NA61/SHINE experi- 
ment [ll|, [13 ' which cover most of the kinematic region 
of interest. The predicted integrated muon neutrino flux 
in the chosen fiducial volume for our data exposure is 
2.02x10^2 cm-2. 

For the first published T2K oscillation analyses, the 
uncertainty on the predicted neutrino flux for this beam 
was as large as 20% [isl. [T^. For this work, however, the 
tuning of the particle production model to NA61 /SHINE 
measurements led to a substantial reduction in systematic 
errors in the flux. With the latest results released by the 
NA61/SHINE collaboration on the kaon production cross 
section based on the 2007 data [l^l, the uncertainty of the 
integrated flux is now about 11%. 

The parameterization of the flux uncertainties is de- 
scribed by normalization parameters in bins of neutrino 
energy and flavor at the near detector. The different 
sources of uncertainty can be separated into two cate- 
gories: the hadron production uncertainty and the T2K 
beamline uncertainty. 




Neutrino Energy (GeV) 

FIG. 2. The ND280 flux prediction with systematic error 
bars, for each neutrino species. The prediction takes into ac- 
count the correct relative fractions of 2010 and 2011 beam 
conditions. 



The uncertainties on hadron production are mainly 
driven by the NA61/SHINE measurements and the 
Eichten and AUaby data [l^ [l^ , and constitute the dom- 
inant component of the flux uncertainty. They include 
the uncertainty in the production cross section, the sec- 
ondary nucleon production, pion production multiplicity 
and kaon production multiplicity. 

The second category of flux uncertainties is associated 
with operational variations in the beamline conditions 
during the data taking. They include uncertainties in 
the proton beam positioning, the off-axis angle, the horn 
absolute current, the horn angular alignment, the horn 
fleld asymmetry, the target alignment, the position de- 
pendence of the flux in the near detector and the proton 
beam intensity. The last two uncertainties were found to 
be very small and are therefore considered negligible. Ta- 
ble |T] shows the contribution of each source of uncertainty 
to the total uncertainty. 



TABLE I. The contribution of each source to the total muon 
neutrino flux uncertainty. 



Error source 


Error (%) 


Production cross section 


6.4 


Secondary nucleon production 


6.9 


Pion multiplicity 


5.0 


Kaon multiplicity 


A O 

0.8 


Off-axis angle 


1.6 


Proton beam 


1.1 


Horn absolute current 


0.9 


Horn angular alignment 


0.5 


Horn field asymmetry 


0.3 


Target alignment 


0.2 


Total 


10.9 



C. The off-axis ND280 detector 

The ND280 detector is a magnetized particle tracking 
device. The active elements are contained inside a large 
magnet, which was previously used for the UAl and NO- 
MAD experiments at CERN. Inside the upstream end of 
this magnet sits a 7r° detector (P0D) consisting of track- 
ing planes of scintillating bars alternating with cither a 
water target and a brass foil, or a lead foil. Downstream 
of the P0D, the tracker consists of three gas-filled time 
projection chambers (TPCs) and two fine-grained de- 
tectors (FGDs) [l^ made up of finely segmented scin- 
tillating bars. The tracker is designed to measure neu- 
trino interactions in the FGDs. The P0D, TPCs, and 
FGDs are all surrounded by an electromagnetic calorime- 
ter (ECal) which detects 7-rays that fail to convert in the 
inner detectors, while the return yoke of the magnet is 
instrumented with scintillator to aid in the identification 
and range determination for muons that exit out the sides 
of the off-axis detector. Fig.[3]shows an exploded view of 
ND280. 

The main detectors used for the analysis presented in 
this paper are the FGDs and TPCs. For this measure- 
ment, the most upstream FGD (FGDl) is used as an ac- 
tive target for the neutrino interactions, while the second 
FGD (FGD2) and the TPCs act as tracking detectors. 

The FGDs are constructed from 9.61 mm x 9.61 mm 
X 1864.3 mm bars of extruded polystyrene scintillator, 
which are oriented along two orthogonal directions per- 
pendicular to the beam. The ND280 coordinate system 
is a right-handed Cartesian coordinate system with the 
z axis running through the central axis of the tracker in 
the beam direction and the y axis in the upward vertical 
direction. 

The FGDl consists of 5,760 scintillator bars, arranged 
in 30 layers of 192 bars each, with the layers oriented 
alternately in the x and y directions. The alternating di- 
rections of the bars allows three dimensional tracking of 
the charged particles. The bars are glued to thin sheets 
of a glass fiber laminate (G-10) and to each other using 
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FIG. 3. An exploded view of the ND280 ofi'-axis detector. 

Plexus MA590 adhesive. The scintillator bars consist of 
polystyrene doped with PPG and POPOP plus a thin 
reflective co-extruded Ti02 coating. A wavelength shift- 
ing fiber, read out by an MPPC, is embedded inside each 
bar. The composition of the FGD is carbon (86.0%), hy- 
drogen (7.4%), oxygen (3.7%), titanium (1.7%), sflicon 
(1.0%) and nitrogen (0.1%), where the percentages rep- 
resent the mass fraction of each element. The number of 
nucleons is given by 

T = NaVfV ■ Pscint ■ /aTT (1) 



a=C,0,H,Ti,Si,N 

= 5.50 X lO^^nucleons 

where Na = 6.022 x 10^'^ mol~^ is Avogadro's number, 
Vpv is the fiducial volume (FV) inside FGDl (see Fig. [4]), 
Pscint = 0.963 g/cm'^ is the density of the scintillator in- 
side the fiducial volume including the glue, coating and 
the air gaps between the modules, a runs over the ele- 
ments present in the scintillator, fa is the mass fraction, 
Aa represents the averaged number of nucleons per nu- 
cleus and Ma is the atomic mass. The ratio of protons to 
neutrons is 53.6:46.4. 

In this analysis, the fiducial volume in the beam direc- 
tion begins at the second x-y module (as encountered by 
the beam), leaving the first x-y module as a veto for in- 
coming particles from upstream interactions. The down- 
stream end of the fiducial volume includes the last layer 
of the scintillator. In the x and y directions, a distance 
equivalent to five bars on each side of the FGD lies out- 
side the fiducial volume. Again, this provides a veto for 
particles generated by interactions outside the FV. 

The second FGD is a water-rich detector consisting of 
seven x-y modules of plastic scintillator alternating with 
six 2.5 cm thick layers of water. It is used as a tracking 
detector in this analysis. 



6 



The two FGDs are sandwiched amongst three TPCs. 
Each TPC consists of an inner box, which is fiUed with an 
Ar:CF4:iC4Hio (95:3:2 by mass) drift gas mixture, within 
an outer box containing CO2 acting as an insulating gas. 
The inner (outer) wahs are made from composite panels 
with copper-clad G-10 (aluminum) skins. The inner box 
panels were precisely machined with an 11.5 mm pitch 
copper strip pattern which, in conjunction with a central 
cathode panel, produces a uniform electric drift field in 
the active drift volume of the TPC, roughly aligned with 
the magnetic field provided by the magnet. 

At each side of the TPCs, 12 micromegas (micro- 
mesh gas detector) modules are arranged in two vertical 
columns. The role of the micromegas is to amplify the 
ionization signals created by charged particles passing 
through the TPCs. 

The micromegas anode is segmented into pads of 
7.0 mm X 9.8 mm (vertical x horizontal) allowing 3D 
track reconstruction of charged particles traversing the 
TPC. The X coordinate is obtained by projecting the ar- 
rival times of the pad signals as discussed in section IV. B. 

The TPCs perform three key functions in the near de- 
tector: 3-dimensional reconstruction of charged particles 
crossing the detector; measurement of the momenta of 
charged particles via curvature in a 0.2 T magnetic field; 
and particle identification (PID) for charged particles via 
ionization (dE/dx). These three functions allow the selec- 
tion of high purity samples of different types of neutrino 
interactions. 



III. SIMULATION OF NEUTRINO 
INTERACTIONS 



of this section. We first describe, in Sec. IIII Al the vari- 
ous interaction models simulated in the event generators. 
The modeling is important to understand the selection 
efficiency and backgrounds, as well as how reconstructed 
and true quantities are related. In Sec. IIII Bl we present 
the systematic uncertainties associated with each kind of 
interaction. In Sec.|Vl these uncertainties arc propagated 
to the reported cross section. 

A. Neutrino interaction model 

In the following, the main interactions related to our 
CC inclusive measurement are discussed. Table lU shows 
the relative fractions of each neutrino interactions type 
simulated in our mixed target. The Monte Carlo predicts 
that 89% of interactions occur on carbon. Some specific 
interactions such as NC elastic scattering, charm pro- 
duction and neutrino-electron clastic scattering are not 
discussed in this paper due to their irrelevance to this 
analysis. We refer to [1, Q for the details. 

TABLE II. Fraction of neutrino interaction processes simu- 
lated in the total FGDl volume by NEUT and GENIE. 



Channel 


NEUT (%) GENIE (%) 


CC quasi-elastic 


39.4 


37.0 


CC Itt resonant production 


18.9 


19.9 


CC coherent n 


1.6 


0.5 


CC other 


11.8 


13.8 


NC 


28.4 


28.7 



Neutrino interactions in the entire detector are sim- 
ulated with the NEUT [| and GENIE H program li- 
braries. While NEUT is used as the primary generator, 
including for the propagation of systematic errors, GE- 
NIE is used for cross-checks and for the generation of fake 
data sets. NEUT and GENIE cover a similar range of 
neutrino energies from several tens of MeV to hundreds 
of TeV and are able to simulate all the nuclear targets 
present in the near detector. 

The neutrino interactions are simulated in both neutral 
current (NC) and CC modes and include the processes 
of elastic scattering, quasi-elastic scattering, single me- 
son production, single gamma production, coherent pion 
production and non-resonant inelastic scattering. Inter- 
actions of hadrons produced inside the nuclear medium 
are also simulated. 

Simulations are a critical aspect of this analysis. In 
the simulated data, neutrino interactions are generated 
outside and within the full detector volume including all 
active and inactive material. This provides information 
necessary to understand the signal as well as backgrounds 
from interactions outside the fiducial volume. The details 
of the simulation process are presented in the remainder 



1. Quasi-elastic scattering 

Both generators use the Llewellyn Smith formalism [l^ 
for quasi-elastic (QE) scattering. In this model the 
hadronic weak current is expressed in terms of the most 
general Lorentz-invariant form factors which do not vi- 
olate G-parity: namely two vector form factors, one 
pseudo-scalar form factor and one axial form factor. The 
vector form factors are measured over a broad range of 
kinematics in electron elastic scattering experiments and 
are fixed by the conserved vector current (CVC) hypoth- 
esis. The pseudo-scalar form factor is assumed to have 
the form suggested by the partially conserved axial cur- 
rent (PCAC) hypothesis [I3. This leaves the axial form 
factor Fa(Q^) as the sole remaining unknown quantity. 
The value of the axial form factor at = 0, Fa{Q), 
is well known from measurements of neutron beta decay 
and the dependence of this form factor can only be 
determined in neutrino experiments. In both generators, 
a dipole form factor is assumed with an axial mass, 
of 0.99 GeV/c2 for GENIE and 1.21 GeV/c^ for NEUT. 

NEUT and GENIE use slightly different approaches in 
the treatment of nuclear effects. Both NEUT and GENIE 
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use the relativistic Fermi gas model (RFG) to describe the 
nuclear environment. The RFG assumes that nucleons 
are bound in a potential well of binding energy, Eb, and 
that they are quasi-free with a maximum momentum set 
to the momentum at the Fermi sphere, pp- The values 
of pp and Eb are measured from electron scattering for 
each target nucleus species individually. While NEUT 
uses the analytical expressions from the Smith- Moniz poj 
model, GENIE incorporates short range nucleon-nucleon 
correlations in the RFG model and handles kinematics 
for off-shell scattering according to the model of Bodek 
and Ritchie f2ll. 



2. Meson production via baryon resonances 



NEUT and GENIE use the Rein-Sehgal [221 model to 
simulate neutrino-induced single pion production. Of the 
18 resonances in the original Rein-Sehgal paper, the 16 
listed as unambiguous in the latest PDG baryon tables 
are included in GENIE [Hi, while NEUT considers all 
18 resonances below the hadronic invariant mass of 2 
GeV/c^. Axial mass values, M^^^ , of 1.21 GeV/c^ and 
1.12 GeV/c2 are used by NEUT and GENIE, respectively. 

For GENIE, the Rein-Sehgal model is used up to a 
hadronic invariant mass of 1.7 GeV/c^, while the analo- 
gous cutoff is 2 GeV/c^ for NEUT. Below this cut value, 
GENIE and NEUT use different descriptions of the non- 
resonant inelastic scattering. This is discussed further in 
Sec. lIIIA4l 



the low parton distributions. The cross sections are 
computed at a fully partonic level and all relevant sea and 
valence quarks are considered. The longitudinal structure 
function is taken into account using the Whitlow param- 
eterizati qn 13011 . The default parameter values are those 
given in [23], which are determined based on the GRV98 
(Glueck-Rey-Vogt-1998) parton distributions. The 
same model can be extended to low energies; and it is 
used for the non-resonant processes that compete with 
resonances in the few-GeV region. 

NEUT and GENIE use different methods to treat the 
non-resonant processes at low W (i.e., low hadronic sys- 
tem invariant mass). Below 2 GeV/c^, NEUT uses a 
W^-dependent function to determine the pion multiplic- 
ity in each interaction. The mean multiplicity of charged 
pions is estimated from data collected in the Fermilab 
15-foot bubble chamber experiment [s^]. In both gener- 
ators, Koba-Nielsen-Olesen (KNO) scaling [s^ is used to 
determine the charged hadron multiplicity. 

Below 1.7 GeV/c^, GENIE uses the Andreopoulos- 
Gallagher-Kehayias-Yang (AGKY) hadronization 
model [sj] to decompose the Bodek and Yang model 
into single pion and two pion production contributions. 
A fraction of these CC-Itt and CC-27r contributions 
are added to the Rein-Sehgal resonance model. The 
fractions are derived from fits to CC inclusive, CC-Itt 
and CC-27r bubble chamber data. The corresponding 
fractions for NC are worked out from the CC fractions 
using isospin arguments. 



3. Coherent pion production 

Coherent scattering results in the production of for- 
ward going pions in both the CC and NC channels. Al- 
though both the NEUT and GENIE Monte Carlos use 
the Rein-Sehgal model [13] for coherent pion production, 
they differ by a factor of 2 in total cross section, as GE- 
NIE uses a recent revision of the Rein-Sehgal model (23| . 
Both are, in turn , sig nificantly larger than recent theoret- 
ical calculations [23. The CC coherent pion production 
of pions has been observed for neutrino energies ranging 
from 2 to 80 GeV [2^. At energies below 2 GeV, obser- 
vations from K2K and SciBooNE are consistent with no 
coherent production [13, . 



4- Non-resonant inelastic scattering 

There are simalirities in how NEUT and GENIE treat 
non-rcsonant inelastic scattering. Both generators use 
the modifications suggested by Bodek and Yang [2^ to 
describe scattering at low Q^. In this model, higher twist 
and target mass corrections are accounted for through 
the use of a new scaling variable and modifications to 



5. Final state interactions 

Final state interactions (FSI) are strong interactions 
affecting particles produced in neutrino interactions as 
those particles traverse the target nucleus. These pro- 
cesses can produce a different observable state as com- 
pared to the state at the initial neutrino-nucleon inter- 
action vertex. In the case of pions, four processes are 
simulated as FSI: QE scattering, absorption, charge ex- 
change and particle production. 

NEUT and GENIE use different microscopic cas- 
cade models to propagate the pion through the nuclear 
medium. In NEUT, the calculation of the probability of 
interaction is separated into low and high pion momen- 
tum regions (< 500 MeV/c and > 500 MeV/c) where 
different models are adopted. In the region < 400 
MeV/c, the A-hole model [11] is employed. For > 500 
MeV/c, np scattering cross sections are used to calcu- 
late the interaction probabilities. For > 400 MeV/c 
and < 500 MeV/c, a linear blending is done to alle- 
viate discontinuities at 500 MeV/c between the A-hole 
model and the scattering data. GENIE makes use of the 
INTRANUKE model [sg, while NEUT uses the model 
described in [3]. The FSI model in both generators are 
tuned to external pion-nucleon scattering data [^, [s^] . 
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B. Neutrino interaction model uncertainties 

In this section, we explain how the cross section un- 
certainties have been calculated. We use a data-driven 
method where the NEUT predictions are compared to 
available neutrino-nucleus data in the energy region rel- 
evant for T2K. We fit the free parameters of the models 
implemented in NEUT, and introduce ad hoc parame- 
ters, often with large uncertainties, to take account of 
remaining discrepancies between NEUT and the data. 



1. Estimation of the CCQE scattering uncertainty 



TABLE III. Summary of the CCQE cross section uncertainties 
used in this analysis, s^^, x^^ and x^^ denote the CCQE 
normalization for different energy ranges, while xsf describes 
the nuclear model applied, where xsf = Q corresponds to the 
use of the relativistic Fermi gas model and xsf = 1 to the use 
of the spectral function. 



Parameter Energy range (GeV) Nominal value Error 





all Ey 


1.21 GeV 


37% 


Xi 


0.0 <Ey < 1.5 


1 


11% 




1.5 <Ey < 3.5 


1 


30% 


J/3 


3.5 < E„ 


1 


30% 


XSF 


all Ey 





100% 


Pf 


all Ey 


217 McV/c 


14% 



The uncertainty in the quasi-elastic cross section is es- 
timated by comparing data from the MiniBooNE exper- 
iment to NEUT. In this comparison, NEUT CCQE in- 
teractions are simulated using the predicted MiniBooNE 
flux 38| and tuned to the MiniBooNE double-differential 
muon CCQE data [H to fit for the best value of M^^ 
and CCQE normalization. The error on each NEUT pa- 
rameter is determined as the difference between the fitted 
value of the parameter and the nominal. For CCQE, we 
take the MiniBooNE data below 1.5 GeV, and assign the 
uncertainty of 11% as reported by the MiniBooNE collab- 
oration [39[. To allow for the discrepancy in CCQE cross 
section between NOMAD ^ at --10 GcV and Mini- 
BooNE at '--I GeV a 30% error has been set above 1.5 
GeV. 

To explore the dependence of the result on the use of 
the simple RFC nuclear model, a spectral function (SF) 
model [41| implemented in the NuWro generator [42| is 
used for comparison. The major contribution to the SF 
comes from the shell model and the remaining -^20% from 
correlated pairs of nuclcons. The latter part accounts for 
a high momentum tail in the nuclcon momentum distri- 
bution, which extends far beyond the Fermi momentum. 
The effective binding energy (equivalent to Eb in the 
Fermi gas model) is on average larger which makes the 
cross section smaller for the SF implementation as com- 
pared to RFC. The fractional difference in CCQE event 
yields calculated with the REG model and the SF is used 
to set the uncertainty in each bin of reconstructed muon 
momentum and angle. Because the uncertainty on the 
Fermi momentum is large enough to have a non- negligible 
effect on the shape of the spectrum for CCQE events, 
a systematic error on the Fermi momentum is added over 
and above the systematic uncertainty assigned using the 
spectral function comparison. In this case, the uncer- 
tainty on the Fermi momentum is taken from electron 
scattering data [2^ and increased to cover MiniBooNE 
data at low Q^. 

Table IIIII summarizes the uncertainties related to the 
charged current quasi-elastic interactions for the NEUT 
generator. 



2. Estimation of the CC non-QE scattering uncertainty 

To constrain the single pion production, we perform a 
joint fit to the MiniBooNE data sets for CC-l7r° [H, CC- 
l7r+ [i^l and NC-l7r° production [i^, since these sets are 
connected in the underlying (Rein-Sehgal) model. Wc fit 
to the reconstructed distributions in the CC channels 
and the pion momentum distribution in the NC channel. 
Nine parameters, described in subsequent paragraphs, are 
included in the fit: M^^^, WcS, the CC coherent nor- 
malization {x'^''^''^"''"), CC-Itt normalizations (xf*^^'^ and 
x^'"^'^), CC-other shape {xccother), NC-Itt" normahza- 
tion (x^'^i'""), NC-l7r=t normalization (x^'^^''*), NC co- 
herent normalization {x^'-^'^"'^'), and NC-other normaliza- 
tion (x^C'o*'''^''). The MiniBooNE data cannot constrain 
all the parameters separately. The fit can, however, re- 
duce the total parameter space through correlations. Be- 
cause the MiniBooNE data sets used in these fits can 
not constrain the parameters xcc other, x^'~^^'^ , x^'-^'^"''- 
and ^ we add penalty terms to prevent large ex- 

cursions away from the generator defaults in the best fit 
values of these parameters. These terms represent con- 
servative prior uncertainty estimates on these parameters 
and are derived with the use of K2K and SciBooNE mea- 
surements. 

The MiniBooNE data directly constrain M^^^ , 
^NCin g^j-^^ CC-Itt normalization parameter for neu- 
trino energies below 2.5 GcV, Xi'~^^'^ . For energies above 
2.5 GeV, we assign a conservative 40% error to the nor- 
malization of CC-Itt production, motivated primarily by 
NOMAD data The Wes parameter modifies the 

width of the hadronic resonance, but not its normaliza- 
tion. It allows an adjustment of the shape of the jp^o | 
spectrum of the NC-Itt*^ channel to improve agreement 
with data. The error on this parameter is taken to be 
50%. 

A 100% error has been set for the CC coherent pion 
production, x'~^'~^''°'^ . This is driven by the fact that K2K 
and SciBooNE p7l.[28| data indicate there is much less co- 
herent charged pion production by neutrinos with energy 
below 2 GeV than predicted by the original models. 
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The xccother parameter modifies a combination of 
otlier CC cross section cliannels as a function of E^. The 
interactions contributing to this category are the CC-utt 
production, which are interactions with more than one 
pion in the final state that have a hadronic mass between 
1.3 GeV and 2 GeV, and DIS or CC resonant interactions 
with rj I Kj 7 production. From external data sets [lO] , the 
uncertainty is known to be of the order of 10% at 4 GeV. 
Using this as a reference point, the error is defined as 
decreasing with the neutrino energy (0.4 GeY / E^, ). 

The NC-Itt^ and the NC-other normalization uncer- 
tainties are set to 30% following the studies done for the 
first published T2K oscillation analyses [l^. In Mini- 
BooNE, there are very few events corresponding to these 
channels and the normalization of these events is not well- 
constrained in the fit. 

A 30% uncertainty on the NC coherent normalization 
factor is used. This conservative estimate is motivated by 
the observation of a 15% discrepancy between the NEUT 
prediction and the SciBooNE measurement of NC coher- 
ent production [i^, together with a 20% systematic error 
in those data. 

In addition to the nine parameters mentioned earlier, 
two additional parameters are considered: the xit^e,, pa- 
rameter and the rate of 7r-less A decay, i^-jess. The 
xittE^ parameter is an empirical parameter that covers 
the discrepancy between the MiniBooNE measurement 
of the CC-l7r+ cross section versus E^ and the NEUT 
prediction using the best fit parameters from the fit to 
MiniBooNE data described above. The discrepancy is as 
large as 50% at E^ = 600 MeV. 

TT-less A decay, also called A re-absorption [l^l, is the 
interaction of the A within the target nucleus prior to 
decay, yielding no pion. This process is assumed to occur 
in ~20% of resonant interactions [1^. It is currently im- 
plemented in NEUT independently of energy and target, 
and results in a CCQE-likc event. An absolute error of 
20% is assigned to this process. 

Tabic Hvl summarizes the different uncertainties on the 
non-quasielastic channels for the NEUT MC. 



3. Estimation of the FSI uncertainty 

In theory, the uncertainties on the FSI parameters (ab- 
sorption, charge exchange, QE scattering and inelastic 
scattering) are correlated with the other cross section pa- 
rameters. In this analysis, however, we assume them to 
be independent, as a first approximation. Therefore, the 
uncertainty on the FSI contribution is added in quadra- 
ture to the other sources in the reported cross section 
(see Sec. IV A|) . The uncertainties on the parameters that 
scale the microscopic interaction probabilities are shown 
in Table |Vl They have been estimated from comparison 
to external tt — C scattering data [SUI- 



TABLE IV. Summary of the CC non-QE cross-section uncer- 
tainties. 



Parameter 


Energy range 


Nominal 


Relative 




[KjG V ) 


value 


Error 




nil P 

all Hii/ 


J LU KjKi V 


y /o 


yy eft 


all Ev 


\ 


/o 


^CCcoh. 


all El, 


1 


100% 






1.63 


26% 


CClTT 

X2 




1 
± 


/o 


XCCother 


all E„ 





40% at 1 GeV 




all E„ 


1.19 


36% 




all £^ 


1 


30% 


^NCcoh. 


all Ey 


1 


30% 


^NCother 


all Ey 


1 


30% 




all El, 


off 


50% 


'^TT — leSS 


all Ey 


0.2 


100% 



TABLE V. A list of parameters used to calculate the FSI 
uncertainties for the NEUT generator only. The low and high 
momentum range refer to the pion momentum smaller and 



greater than 500 MeV/c (see Sec. IIII A 5|l . 

Parameter Error 

Absorption (low momentum) 45% 

Charge exchange (low momentum) 60% 

QE scattering (low momentum) 60% 

Charge exchange (high momentum) 30% 

QE scattering (high momentum) 40% 



Inelastic scattering (high momentum) 50% 



IV. RECONSTRUCTION AND SELECTION OF 
CHARGED CURRENT NEUTRINO 
INTERACTIONS 

At the core of this analysis is the reconstruction and 
event selection in the MC and the data, as described in 
detail in this section. In this analysis, we select CC inclu- 
sive candidates by identifying negatively charged muon- 
candidate tracks in the tracker. The number of selected 
events in reconstructed bins is used to infer the number 
of true events and, consequently, the CC inclusive cross 
section, as described in Sec. El The selected events in 
each bin are dependent on the reconstruction and detec- 
tion efficiency. Below, in Sec. I IV A) we first review how 
the various pieces of MC simulation are included in the fi- 
nal MC. Sec. II V B] explains in more detail how the tracks 
are reconstructed in the TPCs and FGDs. The event 
selection, described in Sec. IIVC| relies on these recon- 
structed tracks. The performance of the event selection 
and reconstruction is provided in Sec. lIVDI along with a 
data-MC comparison. Finally, a discussion of the detec- 
tor response systematic uncertainty is given in Sec. IIVE| 
while Sec. |V] describes how the detector uncertainties are 
then propagated to the final cross section. 
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A. Monte Carlo simulation 

The MC simulation of neutrino interactions in the near 
detector can be divided into several steps. First, the neu- 
trino flux is simulated and propagated to the near detec- 
tor (see Sec. Ill Bp . The neutrino interactions in ND280 are 
then simulated by the two generators NEUT and GENIE 
(see Sec. IIII[) . while only NEUT is used to simulate in- 
teractions located outside ND280, e.g. the pit wall, sand 
and magnet. The secondary interactions in ND280, the 
response of the active detector components and readout 
electronics are simulated using a mix of GEANT4 (ioj 
and custom software. The MC statistics used for this 
analysis correspond to 1.7x10^^ POT for interactions in 
ND280 (for both NEUT and GENIE) and T.OxlQi^ POT 
for interactions outside ND280. 



B. Track reconstruction 

The TPC detector was previously discussed in 
Sec. Ill CI The first step in the TPC reconstruction is the 
application of the gain calibration constants and the re- 
moval of dead and noisy pads, resulting in a waveform 
representing the charge acquired in a single pad as a 
function of the readout time. The charge and time of 
pulses observed in the waveform are extracted, and clus- 
ters are formed by joining coincident pulses on adjacent 
pads in vertical columns. Contiguous clusters are then 
joined to form track candidates and the track parameters 
are obtained with a maximum-likelihood fit to the ob- 
served charge distribution assuming a helical trajectory 
and accounting for the diffusion of the electrons ionized 
by the track. The drift distance is calculated by extract- 
ing the time for the passage of the track (to) by matching 
the track to objects reconstructed in the fast scintillator- 
based detectors, i.e., the FGDs, POD and ECals. 

To perform the particle identification (PID) on a track, 
the charge in the clusters is first corrected for the dis- 
tance traveled by the ionization electrons. This distance 
is estimated based on the reconstructed track position. 
A truncated mean charge is then formed from the cor- 
rected cluster charges. The expected charge deposition 
for each of several particle hypotheses is calculated for 
the measured momentum of the track, and compared to 
the measured charge. 

The FGD data is reconstructed after the TPC recon- 
struction. First the FGD hits are separated into various 
clusters in time. The TPC tracks are then matched to 
FGD hits to form a three dimensional track. 



C. Event selection 

The event selection consists of a series of cuts designed 
to select CC inclusive interactions in FGDl. Interac- 



tions in FGD2 are not considered in this paper, as they 
involve a more complex target consisting of a mix of car- 
bon and oxygen. An extension of this analysis to FGD2 
interactions is anticipated in the future. In addition, 
there is no attempt in this analysis to select backward- 
going muons. Therefore, the CC inclusive selection cuts 
that follow are based on the observation of forward-going 
tracks compatible with being negatively charged muons. 
The cuts used in this analysis are described below. 

1. Data quality flag 

We require that the whole ND280 off-axis detector 
in a full spill was working properly. 

2. Time bunching 

Tracks are grouped together in bunches according 
to their times. This treats neutrino interactions in 
two different bunches within the same beam spill as 
two different events, reducing the accidental pile-up 
of events. The bunch width is ~ 15.0 ns in data. 
Tracks are associated within a bunch if they deviate 
from the mean bunch time by less than 60 ns (i.e. 
4 times the bunch width in data), other tracks are 
removed. 

3. Negatively charged track in the FGDl fiducial vol- 
ume 

We require at least one negatively charged track 
(with FGD and TPC components) starting inside 
FGDl's fiducial volume with more than 18 TPC 
clusters. The interaction vertex is defined as the 
beginning of this track described by the coordi- 
nates {xo,yo, zq). It corresponds, in general, to the 
place where the 3d-fittcd track intercepts the verti- 
cal plane of the upstream-most matched FGD hit. 

The FGDl volume and fiducial volume are shown 
in Fig. m As mentioned earlier, five bars on either 
end of each layer in the FGD are excluded from the 
fiducial volume in the x and y dimensions, while the 
upstream z cut places the fiducial volume just after 
the first xy module. The fiducial volume contains 
therefore 14, xy modules in which the x and y layers 
contain 182 scintillator bars. With this definition, 
the fiducial volume cut becomes |a;o| < 874.5 mm, 
\yo - 55 mm| < 874.5 mm, and zq e [136.9,447.0] 
mm. 

The requirement that the track should contain at 
least 18 clusters is called the TPC track quality cut 
and it rejects short tracks for which the momentum 
reconstruction and particle identification is less reli- 
able. The choice of this particular value of the qual- 
ity cut is based on studies of the kinematic bias for 
tracks of different length. Since only a small frac- 
tion of the selected tracks have fewer than 19 hits, 
the effect of the quality cut, and the systematic er- 
ror associated with it, is very small. 
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FIG. 5. Side view of a charged current event candidate in the 
tracker region of the near detector, which shows the bend- 
ing due to the magnetic field. The muon candidate is recon- 
structed with an angle of 40° and a momentum of 566 MeV/c. 



D. Performance 

Table IVII shows the number of selected events in the 
data and the prediction of the two generators, after each 
cut. 



FIG. 4. Schematic view of FGDl and its fiducial volume (FV) 
drawn in red. In the y-coordinate, the detector is displaced 55 
mm upwards relative to the center of the ND280 coordinate 
system, which is centered in the magnet. 

If there is more than one negatively charged track 
passing these cuts, we select the highest momentum 
track as the muon candidate. 



TABLE VI. The number of selected events for the different 
cuts as predicted by the NEUT and GENIE generators com- 
pared to data. The numbers are normalized to the data POT 
of 10.8x10". 



Cut 


Data NEUT GENIE 


Good negative track in FV 
Upstream veto 
PID cut 


8837 8899 8673 
6243 6582 6351 
4485 4724 4536 



4. Upstream veto 

The goal of this cut is to remove events entering 
the FGDl fiducial volume from the upstream face 
of the detector. Events with more than one recon- 
structed track are rejected if there is a track starting 
more than 150 mm upstream of the muon candidate 
starting position. 

5. TP C particle identification (PID) 

Given the estimated momentum of the muon candi- 
date, a discriminator function is calculated for the 
muon, pion, and proton hypotheses based on the en- 
ergy loss of the track in the TPC. This cut, applied 
to the muon candidate, rejects electrons at low mo- 
mentum (below 500 MeV/c) and removes protons 
and pions. 

Taken together, the cuts above define the CC inclu- 
sive selection in FGDl. The events surviving these cuts 
are included in the final data sample for further analysis. 
Fig. [5] shows one of the events selected via these cuts. 



In Table I VIII the resulting inclusive charged current 
efBciency and purity are shown for each cut. The pri- 
mary effect of the upstream veto cut is to remove events 
taking place in the P0D or in magnet coils that cause a 
secondary interaction in the FGD. The increase of purity 
after the PID cut is due to the PID discriminator's ability 
to separate muons from the low momentum electrons pro- 
duced outside the fiducial volume by mostly NC events. 
As a consequence of the selection, the muon candidate in 
the signal sample is identified correctly 96% of the time. 



TABLE VII. The CC efliciency, e, and purity, V, estimated 
with the MC simulation. 



Cut 


NEUT 


GENIE 


£ (%) V (%) 


e (%) V (%) 


Good negative track in FV 


56.8 52.7 


58.3 51.1 


Upstream veto 


54.4 68.4 


56.2 67.4 


PID cut 


49.5 86.8 


51.2 86.2 



In Figs. [S] and [71 the efficiency as a function of 
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muon momentum and angle are shown as estimated from 
NEUT and GENIE MC simulations. The efficiency for 
backward-going muons is very low since the reconstruc- 
tion does not attempt to reconstruct backward-going 
tracks. The non-zero efficiency for backward-going muons 
arises from the reconstruction of events with a forward- 
going pion possessing a negative charge that is misiden- 
tified as a negatively charged muon, while the real muon 
is going backwards. This effect is included in the cross- 
section determination of Sec. IV Bl 



TABLE VIII. The background composition of the CC inclu- 
sive selection according to the NEUT and GENIE MC gener- 
ators. 



Type 


Fraction of sample (%) 
NEUT GENIE 


Outside FV but in FGDl 

Outside FGDl 
Neutral Currents in FV 
CCue in FV 
'^M,(e) in FV 


0.94 ± 0.14 0.91 ± 0.14 

8.16 ± 0.40 8.36 ± 0.40 

3.17 ± 0.26 3.59 ± 0.28 
0.27 ± 0.08 0.26 ± 0.08 
0.68 ± 0.12 0.71 ± 0.12 


Total 


13.2 ± 0.5 13.8 ± 0.5 
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FIG. 6. The event selection efficiency as a function of the 
muon momentum with its statistical error bars. 



o 
c 

1 0.8- 
w 



0.6 
0.4- 
0.2- 




-NEUT 



■ GENIE 



-0.5 



0.5 1 

True cos 9,, 



FIG. 7. The event selection efficiency as a function of 
the muon angle with its statistical error bars. The non- 
zero efficiency for the backward-going muon is due to mis- 
reconstruction. 

The various backgrounds in the selected sample are 
shown in Table I Villi The main background comes from 
interactions outside the FGD (external background). 

The phase space for events selected in the data is shown 
in Fig. [HI The momentum and angular distributions are 
shown in Fig. [HI where the level of the various back- 
grounds is taken from the NEUT Monte Carlo and the 




500 1000 1500 2000 

Reconstructed (MeV/c) 



FIG. 8. The distribution of selected events in the data in the 
(p^,cos6lp) plane. 



E. Detector response uncertainties 

The systematic errors on the detector response can 
be separated into two categories: the uncertainties on 
the backgrounds in the selected sample and the uncer- 
tainties in the reconstruction. Depending on the un- 
certainty source, the systematic error is computed by 
data/MC comparison, the use of special Monte Carlo 
calculations, the study of cosmic ray muons, or exter- 
nal measurements. For each source, a covariance matrix 
on the predicted number of events is computed according 
to the two-dimensional binning in momentum and angle 
described in Sec. lV Al The two main uncertainties due to 
the detector response come from the momentum uncer- 
tainty arising from distortions in the magnetic field and 
from external backgrounds. 

The magnetic field was mapped for different magnet 
currents in 2009 and 2010 [l|. In addition, a pattern of 
thin aluminum discs and strips were glued to the cathode 
in the TPC. By the use of a laser emitting 266 nm light. 
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Reconstructed cos 

FIG. 9. The muon momentum and angle distribution for the 
selected events in the data and MC. The GENIE prediction is 
shown by the solid blue line, while the backgrounds and signal 
as derived from NEUT are shown in filled colors. 



electrons from the aluminum are ejected and measured in 
the TPC [l3l. The displacement of the electrons is used 
to estimate the magnetic field distortions. Depending on 
the bin, the size of this uncertainty varies from 0.3% to 
7% on the predicted number of events. 

The external background uncertainty is evaluated by 
dividing it into two main categories: the external back- 
ground coming from outside FGDl and the background 
coming from inside FGDl but outside the fiducial vol- 
ume. A 20% uncertainty is assigned based on a com- 
parison of the interaction rate of the background outside 
FGDl betvifeen data and Monte Carlo. A reconstruc- 
tion uncertainty is also assigned for several categories of 
events such as those containing tracks with very high an- 
gles, tracks in which two FGD hits are missing or tracks 
where the matching reconstruction has failed. The sys- 
tematic uncertainty is assigned based on the difference in 
the failure rate between the data and MC for each cate- 
gory of reconstruction considered. The total uncertainty 



assigned to the external background varies from 0.4% to 
9%, depending on the bin. 

A momentum scale uncertainty, determined using the 
aforementioned magnetic field measurements, is included 
in the analysis. The momentum resolution uncertainty 
is estimated by comparing the momentum difference in 
different TPCs in data and MC for tracks crossing at least 
two TPCs. 

Reconstruction uncertainties associated with TPC- 
FGD matching efficiency, tracking efficiency and hit effi- 
ciency are also evaluated. The TPC-FGD matching effi- 
ciency is calculated for each bin separately using control 
samples such as cosmic ray tracks passing through TPC2 
and FGDl. The error on this efficiency varies from 0.2% 
to 2%, leading to an uncertainty on the predicted number 
of events smaller than 1% in all bins. The uncertainty on 
the tracking efficiency is determined by comparing tracks 
with high muon purity crossing the entire detector in data 
and MC. It is found to be 0.5%, leading to an uncertainty 
smaller than 0.8% on the predicted number of events. 
The hit efficiency uncertainty is estimated to be 0.1% by 
comparing the distribution of the number of clusters for 
data and MC. This uncertainty leads to one of the smaller 
uncertainties on the predicted number of events (below 
0.002%). 

An additional source of error arises when the charge of 
a track is not determined properly. A charge confusion 
probability is computed for data and MC using a sample 
of events with tracks starting in the P0D and crossing at 
least two TPCs by looking at the fraction of events for 
which the charges differ in the different TPCs. The 0.3% 
data/MC difference leads to less than 1.1% uncertainty 
on the predicted number of events in each bin. 

The uncertainty on the particle identification is ob- 
tained by using samples with a high purity of muons. 
Generally, these are tracks starting in the P0D which 
cross the three TPCs. For such tracks, the difference of 
the energy loss from the expectation is computed for data 
and MC. The uncertainty is then based on the difference 
between data and MC leading to an uncertainty smaller 
than 0.6% in all bins. 

Finally, the uncertainties due to pile-up events arising 
from neutrino interactions in the near detector and in the 
sand are estimated to be 0.2% and 1.5%, respectively. 
The uncertainty from cosmic ray-induced pile-up events 
is estimated to be 0.1%. The effect of cosmic rays on the 
number of predicted events is negligible. 

Table IIXI shows all the detector uncertainties studied 
for this analysis. 



V. FLUX-AVERAGED CHARGED CURRENT 
CROSS SECTION ANALYSIS 

In what follows, the methodology for extracting the 
cross section from the selected events is described. 
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TABLE IX. A summary of all the systematic errors associated 
with the reconstruction and backgrounds. The first column 
lists all the uncertainties taken into account for the CC anal- 
ysis, the second column presents the samples used to estimate 
the uncertainty and the last column shows the error size on 
the predicted events, depending on the bin. The uncertainties 
due the various backgrounds are shown at the bottom of the 
table. 



Systematic Error 



Data Sample Error (%) 



TPC momentum distortion 

TPC momentum scale 
TPC momentum resolution 
TPC-FGD matching eff. 
TPC track eff. 
Hit efficiency 
Charge mis-ID 
TPC Particle ID (PID) 



special MC 0.3 - 7 

external data 0.1 — 2.4 

Beam data/MC 0.2 - 2.3 
sand muon -f- cosmics 0.2 — 1 

Beam data/MC 0.05 - 0.8 

Beam data/MC < 0.002 

Beam data/MC 0.2 - 1.1 

Beam data/MC 0.02 - 0.6 



External background 
Sand muon background 
ND280 pile-up background 
Cosmic ray background 



several samples 
special MC 

Beam data/MC 
special MC 



0.4-9 
0.1 - 1.1 

0.2 
negligible 



TABLE X. Reconstructed muon momentum and angle bin- 
ning with the correspondence to the one-dimensional binning 
index. 



cos 


Reconstructed index number 


[0.94, 1] 


3 


7 11 15 


19 


[0.9,0.94] 


2 


6 10 14 


18 


[0.84,0.9] 


1 


5 9 13 


17 


[-1,0.84] 





4 8 12 


16 



P^(GeV/c) 



[0,0.4] [0.4,0.5] [0.5,0.7] [0.7,0.9] [0.9,30] 



TABLE XI. The true muon momentum and angle binning 
with the correspondence to the one dimensional binning index. 



cos 6^ 


True index number 


[0.94, 1] 


4 


9 


14 


19 24 


[0.9,0.94] 


3 


8 


13 


18 23 


[0.84, 0.9] 


2 


7 


12 


17 22 


[0, 0.84] 


1 


6 


11 


16 21 


[-1,0] 





5 


10 


15 20 


p^(GeV/c) 


[0,0.4] [0.4,0.5] [0.5,0.7] [0.7,0.9] [0.9,30] 



Sec. IV Al summarizes how the central values are deter- 
mined, while Sec. IV Bl provides the details of how the 
uncertainties were determined. 



A. Method 

To calculate the flux-averaged CC cross section, we use 
a method based on Bayes' theorem [H^l to unfold the 
number of reconstructed and selected events in each mo- 
mentum and angle bin. The result of this unfolding gives 
the number of inferred events Nk in "true" bins A:, while 
the number of selected events in a reconstructed bin j 
is Nj'^^. As demonstrated in [H^l, reconstructed and true 
bins do not necessarily have to be the same. As seen 
previously, this analysis does not attempt to reconstruct 
backward-going muons and all events in the cos inter- 
val between —1 and 0.84 have been placed into a single 
bin (see Table |X|. As shown in Fig. [7l there is almost 
no efficiency for the backward-going angle. We therefore 
split the true cos 9 bins into a truly backward bin [-1,0] 
and a more forward bin [0,0.84]. We will determine the 
double differential cross section in the forward direction 
only, while the total cross section is extrapolated into the 
backward direction using the MC and Bayes' theorem. 
The binning and its association to the one-dimensional 
index is shown in Table IxH 

The flux-averaged cross section in the one-dimensional 
bin, k, defined in Table IXIl is given by 



T<f) T<f) 



1 1 



(2) 



where T is the number of target nucleons, is the in- 
tegrated fiux, N^. is the number of inferred events, Ujk 
is the unfolding matrix, Nj°^ is the number of selected 
events in the reconstructed bin j, Bj the number of se- 
lected background events in this bin as predicted by the 
MC simulation. The parameter is the number of re- 
constructed bins and Ck is the efficiency in the true bin k 
predicted by the MC simulations. 

The unfolding matrix gives the probability that an 
event was created in bin k given that it was reconstructed 
in bin j, and is defined according to Bayes' theorem as. 



Ujk = P{k\j) 



PU\k)Poik) 
Y,P{j\a)P,{a) 



(3) 



where nt is the number of true bins and Po{k) is the 
probability to have a CC interaction in the true bin k, 



Po{k) 



No 



(4) 



where TVtot is the total number of events generated and 
Nk the ones generated in the true bin k. P^ik) \s given in 
Table IXIII together with the efficiency in each bin. Note 
that the model dependences entering the analysis are con- 
centrated mainly in the background prediction and in the 
probability Po{k). 

P{j\a) and P{j\k) are the probabilities to have an 
event reconstructed in the bin j when it has been gener- 
ated in the true bin a or fc, respectively. P{j\k) is esti- 
mated using the MC simulation as described in Sec. IfV Al 
and is defined as the number of CC events reconstructed 
in bin j and generated in bin fc, Sjk, divided by the num- 
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TABLE XII. The efficiency, e^, and probability distribution 
for CC events, Pq (fc) as simulated by the nominal NEUT MC 
for each true momentum and angle bins. 



(GeV/c) 


cos 


ek (%) Po(fc) (%) 


[0.0, 0.4] 


[-1,0] 
[0, 0.84] 


1.2 
26.0 


13.9 
18.5 




[0.84, 0.90] 


62.1 


1.1 




[0.90, 0.94] 
[0.94, 1] 


60.3 
56.0 


0.7 
0.7 


[0.4,0.5] 


[-1,0] 
[0, 0.84] 


3.0 
45.6 


0.9 
9.3 




[0.84, 0.90] 


78.1 


0.9 




[0.90, 0.94] 
[0.94, 1] 


83.1 
84.2 


0.5 
0.6 


[0.5,0.7] 


[-1,0] 
[0, 0.84] 


7.2 
55.1 


0.1 
10.5 




[0.84, 0.90] 


78.4 


2.1 




[0.90, 0.94] 
[0.94, 1] 


82.4 
85.5 


1.4 
1.5 


[0.7,0.9] 


[-1,0] 

[U, U.o4J 


28.3 
61.7 


0.0 
3.5 




[0.84,0.90] 


74.2 


1.3 




[0.90, 0.94] 


79.3 


1.0 




[0.94, 1] 


87.5 


1.3 


[0.9,30.0] 


[-1,0] 


0.0 


0.0 




[0, 0.84] 


63.9 


4.0 




[0.84, 0.90] 


73.4 


3.5 




[0.90, 0.94] 
[0.94, 1] 


76.6 
75.0 


4.2 
18.6 



ber of interactions generated in the true bin fc, Nk, 

pim = f . (5) 

Note that contains all the CC events that were cor- 
rectly selected as well as those that were missed by the 
selection. 

Applying the definitions in Eqs. [3l |4]and[5l the unfold- 
ing matrix can be re- written as, 

Ujk - (6) 
and is shown in Fig. [TOl 

B. Cross section uncertainties 

1. Statistical error 

The statistical error is computed with events simulated 
using NEUT. The nominal NEUT MC is fluctuated 1000 
times, following Poisson statistics. The fluctuation is 
done simultaneously for TVJ''', Sjk, Bj and the number 
of missed events needed for the efficiency. In this case, 
iV|°' is the number of events selected as predicted by the 
MC, scaled down to the data POT and fiuctuated ac- 
cordingly. The error is then the RMS of the fluctuated 



Unfolding Matrix for CC selection (NEUT) 




Reconstructed bin 



FIG. 10. The unfolding matrix obtained with the NEUT MC 
simulation. The white line corresponds to backward-going 
muons in the highest energy bin, for which our MC prediction 
is zero events. The x-axis and y-axis bin numbering is defined 
in Table 1X1 and Table IXll respectively. 



sample distribution, and contains the effect of the lim- 
ited MC statistics. The size of the error varies from 4% 
to 11%, depending on the bin, and the corresponding co- 
variance matrix is given by 

^ 1000 

= E (^fc - - ^r"), (7) 

s=0 

where is the cross section result obtained following 
Eq. [5] for the s-th fiuctuated sample, and cr^°™ the result 
obtained with the nominal MC. In addition, this method 
has been checked with analytic calculations. 



2. Systematic error 

The sources of systematic error on the cross section 
are the fiux, neutrino interaction modeling, final state in- 
teractions, detector response, unfolding method and the 
knowledge of the number of target nucleons. The model- 
ing of the final state interactions is treated here as inde- 
pendent of the rest of the neutrino interaction modeling, 
as discussed in Sec. IIII B 31 

The propagation of the systematic error of the first four 
sources listed above is done by reweighting the NEUT 
MC. The correlations inside each source of systematic 
uncertainties are taken into account by generating a cor- 
related set of systematic parameters that are used to 
rewcight the MC. This procedure is repeated with 200 
different sets of parameters (throws) and the RMS of the 
difference between the cross-section result obtained for 
the nominal MC and the 200 throws is used to define a 
covariance matrix. The corresponding covariance matrix 
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is given by 



Vki 



^ 200 

s=0 



(8) 



which is different for each source of systematic error and 
(7^; is the cross section result obtained fohowing Eq. [5] 
for the s-th throw. The fractional covariance matrices, 
Vfc;/(cr^°'"cr"°"'), due to the flux uncertainty and the rest 
of the sources arc shown in Fig. 1111 



(|) covariance matrix 



« 25 
c 

S 20 

H 



15 
10 



10.02 
,0.015 



-0.01 



0.005 



5 10 15 

All except flux covariance matrix 



20 25 
True bins 




20 25 
True bins 



where (crfc) is the mean of the distribution in bin k. The 
error due to the unfolding is generally below 1% except 
for the backward bins. 

The last uncertainty that needs to be taken into ac- 
count is the one on the number of target micleons, T, as 
explained in detail in Eq. [1] of Sec. Ill CI The resulting 
uncertainty of 0.67% is then added in quadrature in all 
bins. Table IXIIII summarizes the uncertainties for each 
bin. 

The uncertainty on the cross section due to the flux is 
similar in each bin, and is dominated by the integrated 
flux normalization uncertainty of 10.9%. 

The errors due to the neutrino interaction modeling are 
very different from one bin to another. Backward-going 
events contain the largest systematic uncertainty since we 
have a very low efficiency for this part of the phase space 
and the result is effectively a model-dependent extrapola- 
tion from the forward direction. The neutrino interaction 
modeling error is dominated by the uncertainty assigned 
to the spectral function corresponding to most of the sys- 
tematic error at momenta between 700 MeV/c and 900 
MeV/c and cos 6*^ > 0.84. This behavior is effectively 
expected as demonstrated in [4lj. For the other bins, the 
uncertainty on is also important and dominant for 

momenta below 400 MeV/c. 

The systematic error on the cross section due to the 
uncertainty on the detector response is bigger at low mo- 
mentum, and is dominated by the external background 
uncertainty. For the other bins, the main contributions 
are shared more or less equally between the uncertainty 
in the magnetic field, FGD-TPC matching, charge con- 
fusion and the fiducial mass. 

The systematic uncertainty due to final state interac- 
tions is relatively small, as its main effect is to change 
the number and energy of the pions, which are occasion- 
ally identified as the muon candidate. These constitute 
a small fraction of the sample and therefore the effect is 
small. 



VI. RESULTS 



FIG. 11. The fractional covariance matrix obtained from the 
uncertainties on the flux (top), uncertainties due to all sources 
except the flux (bottom). The x-axis and y-axis bin number- 
ing is deflned in Table IXB 

The uncertainty due to the unfolding method, the fifth 
source of systematic error listed above, is estimated us- 
ing the same fluctuated sample as for the statistical error. 
The difference between the nominal cross section and the 
average cross-section extracted from these 1000 fluctu- 
ated sets is used as the estimate of the uncertainty due 
to the algorithm bias. The covariance matrix for this 
uncertainty is given by 



(9) 



In this section, the flux-averaged double differential 
cross section is presented. In particular, we report 
the measurement for the forward bins in Sec. IVI Al In 
Sec. IVI Bl we include the backward bin to give the flux- 
averaged total cross-section measurement. 



A. Flux averaged differential cross section 

The flux-averaged double differential cross section is 
calculated as. 



da 



dpndcos 9,1 



l-i- I ap 



T^Ap^^qAcos 6^,^13 



(10) 
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TABLE XIII. A summary of the systematic and statistical errors. The error on the number of target nucleons (0.67%) has 
been added in quadrature to the total systematic error. The "unfold.", "(^", "det.", and "FSI" labels represent the systematic 
uncertainty due to the unfolding method, the beam flux, the detector response and the final state interactions, respectively, 
changed systematically following their respective covariance matrix. The "model" column denotes the infiuence of changing all 
the neutrino interaction modeling parameters and channel rates. 



(GeV/c) 


cos 9^ 


unfold. (%) 0(%) model (%) det. (%) FSI (%) 


syst (%) 


Stat (%) 


tot (%) 


[0.0, 0.4] 


[-1,0] 


0.5 


11.4 


18.0 


2.1 


0.5 


21.4 


2.0 


21.5 




[0, 0.84] 


0.6 


12.8 


5.5 


3.6 


1.2 


14.5 


4.9 


15.3 




[0.84, 0.90] 


0.2 


13.1 


10.8 


2.7 


1.4 


17.3 


9.5 


19.7 




[0.90, 0.94] 


1.0 


14.1 


10.7 


5.0 


3.6 


18.8 


12.3 


22.4 




[0.94, 1] 


0.3 


14.0 


12.9 


4.9 


3.0 


20.0 


14.7 


24.8 


[0.4,0.5] 


[-1,0] 


1.2 


12.0 


39.5 


2.7 


0.9 


41.4 


3.2 


41.5 




[0, 0.84] 


0.2 


11.4 


5.7 


1.3 


0.3 


12.8 


4.2 


13.5 




[0.84, 0.90] 


0.6 


11.4 


5.0 


1.0 


0.4 


12.5 


8.6 


15.2 




[0.90, 0.94] 


0.5 


11.7 


5.4 


1.3 


0.5 


13.0 


10.1 


16.4 




[0.94, 1] 


0.5 


13.1 


7.2 


2.3 


0.9 


15.2 


11.7 


19.2 


[0.5,0.7] 


[-1,0] 


1.2 


12.6 


46.1 


1.9 


0.4 


47.9 


9.0 


48.7 




[0, 0.84] 


0.3 


11.1 


3.8 


1.1 


0.4 


11.8 


3.8 


12.4 




[0.84, 0.90] 


0.3 


10.8 


3.4 


0.8 


0.3 


11.4 


6.2 


13.0 




[0.90, 0.94] 


0.4 


11.0 


5.7 


0.8 


0.4 


12.5 


7.3 


14.4 




[U.y4, ij 


0.0 


11.6 


11.4 


1.1 


0.3 


ID. 4 


( .y 


io.z 


[0.7,0.9] 


[-1,0] 


1.8 


13.5 


148.3 


2.0 


0.6 


149.0 


31.5 


152.3 




[0, 0.84] 


0.4 


11.4 


3.2 


1.1 


0.4 


11.9 


5.1 


12.9 




[0.84, 0.90] 


0.4 


10.9 


5.9 


0.8 


0.2 


12.5 


6.8 


14.2 




[0.90, 0.94] 


0.2 


10.7 


11.1 


1.1 


0.5 


15.5 


7.7 


17.3 




[0.94, 1] 


0.2 


11.0 


17.6 


0.9 


0.4 


20.8 


7.0 


21.9 


[0.9, 30.0] 


[-1,0] 




















[0, 0.84] 


0.2 


11.9 


5.6 


1.4 


0.6 


13.3 


5.4 


14.3 




[0.84, 0.90] 


0.2 


11.3 


2.5 


0.9 


0.3 


11.7 


5.8 


13.1 




[0.90, 0.94] 


0.4 


11.1 


2.3 


0.7 


0.4 


11.4 


5.2 


12.5 




[0.94, 1] 


0.2 


10.9 


2.3 


0.8 


0.3 


11.2 


2.9 


11.6 



where Nap is the number of inferred events for the two- 
dimensional binning in momentum and angle labeled by 
a and /?, respectively, while Ap^^a and A cos 6*^^/3 give 
the respective bin widths. Table IXIVI gives the number 
of reconstructed events in data and the background pre- 
diction needed to calculate the number of inferred events 
as defined in Eq. [2l The number of simulated and pre- 
dicted CC events in the FGD volume using the NEUT 
generator are also shown together with the result of the 
unfolding. 



The measured differential cross section is shown in Ta- 
ble IXVI and Fig. [12] only for the forward bins since the 
results obtained for the backward bins are completely 
model-dependent. In the case of the backward bins, 
we rely entirely on the MC to determine the fraction of 
the selected forward-going events which are due to true 
backward-going events. 



B. Flux-averaged total cross section 

The flux-averaged total cross section is calculated by 
taking all the bins into account, including the backward 
bins, although this analysis has very low efficiency in the 
backward bins resulting in larger statistical errors and 
model dependence. 

The fiux-averaged total cross section is calculated as. 



Tcb 



(11) 



where iV*°* is given in Table IXIVI We obtain 

(o-cc)^ = (6.91 ± 0.13{stat) ± OM{syst)) 

(12) 

nucleon 

which agrees well with the MC predicted values from 
NEUT and GENIE that are, 



The cross section was also calculated unfolding with 
the GENIE generator. The results are consistent within 
the errors estimated from modeling uncertainties. 



= 7.27x10- 



-39 



ag™iE)^ = 6.54 X 10 



-39 



nucleon 

cm^ 
nucleon 



(13) 
(14) 
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FIG. 12. The CC-inclusive differential cross section in cm^/nucleon/MeV, with statistical and systematic errors. Each graph 
corresponds to a bin in muon angle. 



The total cross section result for T2K is shown in Fig. [13] 
together with CC inclusive measurements from other ex- 
periments. 



C. Future improvements 

There is on-going work in T2K expected to reduce the 
model dependence in the neutrino event generator and 
to reduce the flux errors. These studies should lead to 
improved CC inclusive cross-section results in the future. 



VII. SUMMARY 

We have selected a sample of CC inclusive inter- 
actions in the tracker of the T2K off-axis near detector. 
We used a largely model-independent method to extract 
the flux-averaged differential ly^ CC cross section in muon 
momentum and angle. These results represent the first 
such measurement done for CC inclusive interactions on 
carbon at a mean neutrino energy of 0.85 GeV and are 



presented as a two-dimensional differential cross section. 
In addition, these results have been used to calculate a 
flux-averaged total CC inclusive cross section that is con- 
sistent with previous measurements in this energy range. 
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TABLE XVI. The total fractional covariance matrix including statistical and systematic errors. The values are given in (%). 
The bin numbering is defined in Table IXll 
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